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^ ' Abstract 

> : 

Starting from the QCD Schrodinger functional (SF), we define a family of renormalization 
schemes for two four-quark operators, which are, in the chiral limit, protected against mix- 
ing with other operators. With the appropriate flavour assignments these operators can be 
interpreted as part of either the AF = 1 or AF = 2 effective weak Hamiltonians. In view 
of lattice QCD with Wilson-type quarks, we focus on the parity odd components of the op- 
erators, since these are multiplicatively renormalized both on the lattice and in continuum 
schemes. We consider 9 different SF schemes and relate them to commonly used contin- 
uum schemes at one-loop order of perturbation theory. In this way the two-loop anomalous 
dimensions in the SF schemes can be inferred. As a by-product of our calculation we also 
obtain the one-loop cutoff effects in the step-scaling functions of the respective renormaliza- 
tion constants, for both O(a) improved and unimproved Wilson quarks. Our results will be 
needed in a separate study of the non-perturbative scale evolution of these operators. 



1 Introduction 



In the Standard Model, four-quark operators typically arise as effective interaction vertices 
when integrating out the large scale physics associated with the weak interactions. Examples 
are the AB = 2 operators 



which mediate B°-B° mixing or their analogues in the case of kaons. Hence, the quantities 
of interest are matrix elements of four-quark operators between hadronic states, which are 
inherently non-perturbative in nature. On the other hand, the four-quark operators are 
originally obtained in perturbation theory and renormalized at a large scale, using e.g. the 
minimal scheme (MS) of dimensional regularization. If lattice QCD is used for the calcu- 
lation of the hadronic matrix elements, the matching to the perturbative renormalization 
schemes poses a problem, since the scale differences involved are potentially large. A general 
strategy to solve this problem has been proposed some time ago [J: its starting point is 
the definition of an intermediate scheme, where the finite space-time volume is used to set 
the renormalization scale. Finite-size-scaling techniques then allow to step up the energy 
scale recursively until the perturbative regime is reached, where the continuum schemes can 
be safely matched in perturbation theory. Previous applications include the running cou- 
pling 2| an d quark mass |3)4)5j . moments of structure functions jEHUl) and the static- light 
axial current jlUlllj . 

This paper is part of a project to apply this strategy to two particularly important 
four-quark operators with phenomenological applications to B°—B° or K°-K° mixing and 
non-leptonic kaon decays |12H15| . We use the QCD Schrodinger functional to define a family 
of finite volume renormalization schemes, and report our results for the perturbative match- 
ing to more commonly used renormalization schemes. The matching procedure is done in 
two steps: first the Schrodinger functional is regularized on the lattice with Wilson-type 
quarks, and the renormalized operators are related to the standard lattice scheme, defined 
by minimal subtraction of logarithms. Then, using results from the literature, this lattice 
scheme can be related to a continuum scheme, such as dimensional reduction (DRED) or one 
of the minimal subtraction schemes (MS) in dimensional regularization. Since the two-loop 
anomalous dimensions are known in these schemes, the one-loop matching then allows the 
two-loop anomalous dimensions to be inferred in the SF schemes, too. 

The paper is organized as follows: in section 2 we start with the definition of the 
operators and their correlation functions in the Schrodinger functional, which are then used 
to formulate the renormalization conditions. After a short review of perturbative matching 
equations between different renormalization schemes (section 3), we collect the equations for 
the reference schemes (section 4) and report the results of our one-loop computation (section 
5). In view of the corresponding non-perturbative computation with Wilson-type quarks, 
we then discuss perturbative lattice artefacts in the step-scaling functions (section 6) and 
present our conclusions. 
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2 Definitions and setup 



The four-quark operators we would like to renormalize are of the form 

°ll = o - 75)^2)(^37m( 1 - 75)^4) ± W'lT/X 1 ~ 75)^4) ("037m (! - 75)^2)] , (2.1) 

where the subscript LL refers to the Dirac structure of two left-handed currents, and it is 
understood that colour indices are contracted within the quark bilinears in round brackets. 
In order to make contact with phenomenological applications, one just needs to assign the 
physical quark flavours. For instance, with the identifications 

tpi = ip3 = s, ip2 = ip4 = d, (2.2) 

the operator 0£" L vanishes while the matrix elements of 0^ L appear in the K°-K° mixing 
amplitude. Replacing strange by bottom quarks, the same operator mediates B°-B° mixing. 
If instead one identifies 

•01 = s, i> 2 = d, ip3 = ipi = u, c, (2.3) 

one obtains the AS = 1 operators relevant to the AI = 1/2 rule in non-leptonic kaon decays, 
in a framework where the charm quark remains an active degree of freedom. 

The operators (|2.1|) can be decomposed in parity-even and -odd components: 

^LL = ^(V-A)(V-A) = ^VV+AA ~~ ^VA+AV (^-4) 

where V, A refer to the Dirac structure of vector and axial vector currents, respectively. In 
the following we will work with the parity-odd component of the operators, 

°VA+AV = 2 [(^17/^2X^37^75^4) + (^17^75^2X^37^4) ± (V>2 ^ ^4)] , (2.5) 

rather than staying with the product of two left-handed currents (j2.1|) as induced by the 
Standard Model weak interactions. Note that in regularizations that preserve chiral sym- 
metry, parity-even and parity-odd components are related by chiral symmetry and are thus 
renormalized in the same way. However, the situation changes in lattice regularizations with 
Wilson-type quarks: while the parity-odd operators OyA+AV are rnultiplicatively renormal- 
izable their parity-even partners Oyv+AA snare the lattice symmetries with four other 
four-quark operators, which leads to a complicated operator mixing problem. Although this 
problem can be solved non-perturbatively ^7j, it is a source of additional uncertainties and 
systematic errors. In this paper this problem is circumvented by imposing renormalization 
conditions on the parity-odd operator components. Besides its technical advantage this strat- 
egy makes sense also from a practical point of view: as has been demonstrated in |18j . the 
introduction of non-standard chirally twisted mass terms ("twisted mass QCD") redefines 
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the physical parity symmetry so that the parity-odd operator components can play the role 
of operators with either physical parity. 1 As a result, the hadronic matrix elements of oper- 
ators with even physical parity can be obtained from correlation functions which only involve 
the lattice operators 1)2.5)1 . We conclude that our choice to renormalize these parity-odd op- 
erator components is irrelevant for chirally symmetric regularizations, but it is advantageous 
with Wilson-type quarks and does not imply any prejudice on possible phenomenological 
applications. 

2.1 Renormalization conditions 

We now choose the lattice regularization with Wilson quarks, possibly 0(a) improved by 
the Sheikholeslami-Wohlert term in the action [2U]. For unexplained notation we refer the 
reader to ref. |2j. We assume that the bare operators 1)2.5)) are defined locally on the lattice, 
i.e. all quark and antiquark fields are taken at the same space-time point. To formulate 
renormalization conditions for the renormalized operators 

(Or)va+AV = ^VA+AV^VA+AV ( 2 ' 6 ) 

we use the standard set-up of the Schrodinger functional |25|26j as described in [2^. We 
consider generic source fields made up of boundary quarks and antiquarks, 

o 12 [T) = a 6 ^Ci(y)rc 2 (z), (2.7) 
o' 12 [T] = a 6 ^c'i(y)rc 2 (z), (2.8) 

where T is a Dirac matrix that must anticommute with 70 as otherwise the source field 
vanishes. This is due to the projectors P± = |(1 db 70), which are implicit in the boundary 
quark and antiquark fields, 

C(x)=P_C(x), C(x) = C(x)P+, (2.9) 

and similarly for the primed fields. The presence of the projectors limits the possible Dirac 
structures for the quark bilinears at the time boundaries. For instance, it is not possible 
to define scalar quark bilinear sources, unless one introduces a non-vanishing angular mo- 
mentum. However, finite momenta typically increase the lattice artefacts, and lead to poor 
signal-to-noise ratios in numerical simulations, so that we are not going to pursue this further. 

The renormalization conditions will be imposed in the massless theory, so that the 
renormalization constants and the renormalized coupling are quark mass independent by 
construction. In the absence of chirally twisted mass terms, standard parity is an exact 
symmetry of the lattice-regularized Schrodinger functional. In order to renormalize the 

x For an alternative solution using axial Ward identities, see |19j . 



3 



parity-odd operators (|2.5j) we thus need a total source with negative parity which contains 
at least 2 quark bilinear sources. Because of the above mentioned problem with the projectors 
at the time boundaries, we decided to introduce a fifth "spectator quark" and use correlation 
functions of the generic form: 

f^t^Tb^cKxo) = l- 3 (o 21 [r A ]o± A+AV (x)o 45 [r B }o' 53 [rc]) ■ (2.10) 

The corresponding Feynman diagram is given in the left of fig. 1, where the quark lines 
correspond to the boundary-to-volume quark propagators H(x),H'(x) of refs. |22I23| . and 
to the boundary-to-boundary propagator K of ref. [21], which contains an explicit time-like 
link variable from Euclidean time T — a to T. 
We then consider the 5 specific cases 

Ff{x Q ) = F±[ 75 , 75, 75] (so), (2.11) 
1 3 

^0*0) = g Yl E ijkF k bi,'yj,'YkKx ), (2.12) 
i,j,k=l 



1 3 

F t&o) = 3 ^^[75,7^7^(^0), (2-13) 
fc=l 

1 3 

Ft(xo) = oE F± ^'^,7ifc](so), (2.14) 
6 k=i 

1 3 

^(^0) = ^E^bu^slN. (2.15) 

6 k=i 

These correlation functions describe transitions between parity-even and -odd states at the 
Euclidean time boundaries, mediated by the parity-odd four-quark operators. In particular, 
parity-even scalar or axial vector states are produced at the lower time boundary by taking 
appropriate combinations of pseudoscalar and vector states. 

In order to obtain renormalization conditions for the four-quark operators based on 
these correlation functions, we first have to take care of the source field renormalization. As 
the boundary quark and antiquark fields are all renormalized multiplicatively by the same 
renormalization constant [2Z1, this can easily be achieved by forming appropriate ratios of 
correlation functions. More precisely, with the boundary-to-boundary correlators 

fi = -^(Oi2hsP' 21 [ l5 }), (2.16) 
1 3 

fcl = -q^J2(°MOM), (2.17) 

k=l 
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Figure 1: The Feynman diagrams for the four-quark correlation functions F^(xq) and the 
boundary-to-boundary correlators fx, k\ at tree level. Euclidean time goes from left to right. 
The double blob indicates the four-quark operator insertion and the dashed links indicate 
the explicit time-like link variable in the boundary-to-boundary quark propagators. 

we consider the following 9 ratios of correlation functions 

i = l,...,5 (2.18) 

(2.19) 

i = 3,4,5. (2.20) 

In these ratios the renormalization of the boundary fields cancels out, i.e. the renormalized 
ratios are obtained by multiplying with the appropriate four-quark operator renormalization 
constant. A renormalization condition is now obtained by choosing one of the ratios (|2.20|) . 
setting the renormalized quark mass to zero, and specifying T/L, xq/L and the angle 8 of 
the spatial boundary conditions |21| . 

All parameters are then fixed, and L remains the only scale in the system. Then one 
requires 

z va+av ; s(so, an)hf(x ) = hf(x )\ go =o, (2.21) 

where the RHS is the free field theory result. The renormalization constant ^vA+AV-s ^ s thus 
obtained at the renormalization scale // = 1/L, and depends implicitly on all the parameter 
choices made. In practice, we followed refs. |3I4| and always set T = L, xq = L/2 and 9 = 0.5. 
This leaves us with 9 different SF schemes labelled by s = 1, . . . , 9 in eq. (|2.21|) . While there 
are good arguments for choosing T = L |3] , there are in general no a priori criteria for a good 
parameter choice. This can only be judged a posteriori by comparing non-perturbative and 
perturbative data, or by looking at the apparent convergence of the perturbative expansion 
of the anomalous dimensions. As it will turn out, there is considerable variation among our 
9 choices of renormalization schemes. 



hf(xo) 



Ft(xo) 

A 3/2 ' 

Ft(xo) 
fl /2 h 



5 



3 Anomalous operator dimensions in perturbation theory 



Operators and parameters are renormalized at the renormalization scale fj,, which in the SF 
schemes is identified with 1/L. A change of scale is then governed by the renormalization 
group. Let us consider QCD with iVf mass degenerate quark flavours and N colours, and its 
Euclidean correlation functions of gauge-invariant composite operators. Limiting ourselves 
to multiplicatively renormalizable operators, any renormalized n-point functions of such 
operators, 

G n (xi,. ..,x n ) = (Oi(xi) • • • O n (x n )) , (3.1) 
satisfies the Callan-Symanzik equation 

The renormalization group functions f3, r and 7 are the /3-function for the coupling and 
the anomalous dimensions for the quark mass and composite operators respectively. In 
quark mass independent schemes the renormalization group functions only depend on the 
renormalized coupling and have perturbative expansions of the form 

00 

m 9 -° vxv*, (3.3) 

fc=0 

00 

rig) 9 -° -s 2 5><? 2fc , (3.4) 

fc=0 

00 

7(9) 9 -° -9 2 Y,^9 2k . (3.5) 

fc=0 

The coefficients are scheme-dependent in general, except for bo, b\ and do, 70. In the normal- 
ization adopted here we have (see refs. |28I29I3(J] for the coefficients up to k = 3 in the MS 
scheme and for further references): 
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Y N ~l N t \(^r-. (3.0) 



3(iV 2 - 1) (A 2 

do = ^ (4tt) , (3.7) 



h = {|iV 2 -(fiV-iV- 1 )iVf}(4vr)- 



(3.8) 



Denoting the anomalous dimensions for the operators Oy A+AV by their leading order 
coefficients are given by |31I32| 

7o± _ *«£EV)-» ,3.9, 
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The universality of these coefficients can easily be understood by changing to another quark 
mass independent scheme. This amounts to finite renormalizations of the form 



cf = 9yjX g (g), (3.10) 

m' = mX m (g), (3.11) 

0' R = R X (g). (3.12) 

The n-point functions of the primed operators satisfy again a Callan-Symanzik equation of 
the form ()3.2j) . and the respective renormalization group functions are then related as follows 

W) = > ( 3 - 13 ) 

g=g(g') 



d 9 J < 

r'(g') = \r(g) +P(g)^- In X m (g)\ , (3.14) 

I 8 9 ) 9=9(9') 

l'(g') = {i{g)+l3{g)^\nXo{g)\ . (3.15) 
I dg J a=o(gl) 



9=9(9') 

Expanding the renormalization factors in perturbation theory, 

00 

X(g) 9 ~°l + ^X {k) g 2k - (3-16) 

k=l 

one finds that 60,61, and e?0)7o remain indeed unchanged, and for the next-to- leading order 
anomalous dimensions one arrives at 

d[ = d l + 2b xW -doX^, (3.17) 
7i = 7l + 2b X^ -7o^ g (1) - (3.18) 

Hence, if 71 is known from a two-loop calculation in some reference scheme, it can be obtained 
in any other scheme by relating the schemes at one-loop order, thereby avoiding a direct two- 
loop computation. The situation for any multiplicatively renormalizable operator hence is 
the same as with the quark mass, where the two-loop anomalous dimension in the SF scheme 
has been obtained along these lines 0. 



4 Reference schemes 

The two-loop anomalous dimensions 7^ have been computed for a variety of schemes |33H36j . 
The first computation was performed by Altarelli and collaborators using dimensional 
reduction (DRED) |37j. This result was later confirmed by Buras and Weisz |34j . who used 
dimensional regularization with both the naive and the 't Hooft-Veltman definition of 75 
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in D dimensions In the DRED scheme, but with the renormalized coupling defined 

in the MS scheme (this differs from the renormalized coupling used in [S3])) the two- loop 
anomalous dimension takes the form 

*- ~ ^ - » - > ± ( ^ + ^ - ?*) } («) 

For later reference we also quote the corresponding results in the NDR ("dimensional reg- 
ularization with naive 75") and HVDR ("dimensional regularization with 't Hooft-Veltman 
75") schemes, as defined in (31] 

± ATzpi (- ^ | f57 19^ 4 



f iTndr = -^{^^(aF-T^ + ^'JI' 4 '* • (4 - 2) 
7tHVDR - ^{^_ 21 _^ t± (-,^_^)} (4 ^. ( , 3) 

In order to obtain the two-loop anomalous dimensions in the SF schemes, we need the one- 
loop relations between the renormalized operators and coupling constants, 

(^VA+Av)sF = (^VA+Av)dRED "^SF.DREdG?)' (4-4) 
f SF (L) = g\^)X g {g). (4.5) 

Here we have denoted the MS coupling by g, and we assume that the SF coupling has been 
defined for N = 3 colours as in |39I40| , There, also the one-loop coefficient for the matching 
of the couplings has been determined: 

X® = 2b ln(fiL)-^-{c 1> o + c 1>1 N i ), (4.6) 

cl,o = 1.25563(4), (4.7) 
= 0.039863(2). (4.8) 

The one-loop relation between the renormalized operators will be established in two steps, 
by first converting to a lattice renormalization scheme (lat), which is obtained by minimally 
subtracting the logarithms |41j . This yields the one-loop relation 

(OvA+Av) S F = ( C, VA+Av)l a t ;k 'sF,lat(9lat), (4.9) 

which will be discussed in more detail in the next section. On the other hand, the relation 
between the operators in the lat-scheme and DRED has been established in refs. |42|43|44] . 
Defining the finite renormalization factor through 

(°VA+Av)dRED = ( VA+Av)lat' ; t'D R ED,lat(9lat), (4-10) 

the one-loop coefficient is given by 

*DRED,lat = {^1 ± *0 + N^Z^} (^)~\ (4.11) 
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with coefficients 



zi = l - (A 7m + A >75 + 2A El + 1) , (4.12) 
z = - (A 7fl + A 7m75 - Ax - A 75 - 1) , (4.13) 



z-i = -(Ai + A 75 )- (A 7m +A 7m75 + A Ei ). (4.14) 



The A's have been defined in |42l43j and are related to the quark propagator and vertex 
functions of quark bilinears. Though gauge parameter dependent in general, the above linear 
combinations are gauge- independent, and a numerical evaluation yields |45I43I46] . 

. . ([-14.06090(9), 5.7854(1), 8.2755(2)] for c sw = 1, 

[zi, zq, z-i\ = < (4.15) 
[[-17.70704(7), -0.9331(1), 18.6402(2)] for c sw = 0, 

where to this order of perturbation theory c sw = ciw = 0, 1 refers to standard and O(a) 
improved Wilson quarks, respectively. With N = 3, we obtain 

= (-°- 213020(2) (4.16) 
DRED,iat [-0.302956(2) forc sw = 0, 

(1) f-0.286293(2) for c sw = 1, 

DREDilat [-0.291138(2) forc sw = 0. 

Finally, the desired one-loop relation between the SF schemes and the DRED scheme is 
obtained by combining eqs. I|4.9|l . (|4.10j) . which to one-loop order implies 

y±(l) __ „±(1) _ ^±(1) f. 1R X 

SF,DRED _ ^SF.lat DRED,laf K" 1 - 10 ) 

The two- loop anomalous dimensions are then related by formula ()3.18|) . identifying the SF 
and DRED schemes with the primed and unprimed schemes, respectively. 



5 One-loop results 

The perturbative expansion of the finite volume correlation functions is straightforward 
albeit a bit tedious. The technique is well-documented so that we refer to refs. 122124] for 
details concerning the gauge-fixing procedure and the parameter tuning necessary to take 
the continuum limit while keeping the volume fixed in physical units. Here we just describe 
the technical details pertaining to the application at hand. We generated double precision 
data for the one-loop diagrams displayed in figs. 2 and 3, for lattice sizes ranging from 
L/a = 4 to 32, and we took steps of 2 in order to have a lattice coordinate for xq = L/2. We 
generated data for both the O(a) improved (c sw = ciw = 1) and unimproved Wilson quarks 
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Figure 2: Feynman diagrams of the self-energy type. 

(c sw = 0). In the case of the 0(a) improved action, we also included the effect of the 0(a) 
counterterm from the boundary proportional to the improvement coefficient q |21j . We did 
not attempt to improve the four-quark operators, as there are several 0(a) counterterms, 
which renders 0(a) improvement impractical. We note, however, that the local operators 
are 0(a) improved at tree level. Two independent sets of data were generated by two subsets 
of the authors, and perfect agreement up to rounding errors was found. We also checked 
the independence of the gauge parameter, and compared with a numerical simulation with 
gauge group SU(3) at large values of (3 = Q/g^ (e.g. = 80). As disconnected diagrams 
(with respect to the quark lines) start contributing at order g$, the numerical values for 
these diagrams set the scale for the expected accuracy of the comparison. 

Having passed these checks, we obtained the one- loop expressions for the renormalization 
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Figure 3: The remaining Feynman diagrams at one-loop order, which are not of the self- 
energy type. 

constants from the renormalization conditions ()2.21j) . With the notation 

oo 

^VA+AV;s(5o,a/L) = l + ^ 5o 2 "Z±(L/a)W, (5.1) 

n=l 

oo 

Ft(x ) = 5>o 2n if(*o) (n) (5.2) 

n=0 
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(and analogously for fx, k\), we find 




ln/f 



kl + ^ + «4 1 »^ln fc <»' . (5.3) 



k («) t (») •= dm. 



Here, the SF schemes s = 1 to s = 5 correspond to i = 1, . . . , 5 and z = 1, SF scheme s = 6 
translates to i = 2 and z = 0, and schemes s = 7 — 9 are obtained with i = 3, 4, 5 and 
z = 1/3. It is assumed that the one-loop expressions are evaluated at the bare mass tuq = 0, 
i.e. the mass counterterms with (see e.g. |47| ) 

(1) / -0.20255651209 C F (c sw = 1), ^ N 2 - 1 
am}.' = < C F = — , (5.4) 

\ -0.32571411742 C F (c sw = 0), 2A^ 

ensure the condition of vanishing renormalized quark mass to one- loop order. We have also 
added the contribution of the boundary counterterm proportional to c t — 1, as indicated by 
the subscript b. These terms only modify the O(a) cutoff effects and we set them to zero in 
the case of unimproved Wilson quarks. 

In all cases it is expected that the one- loop coefficients Z^(L/a)^ have an asymptotic 
expansion in powers of the lattice spacing of the form 



Z s 



oo 

± {L/a)^ ~ J2 (l )" K + s - HL/a)} • (5.5) 



u=0 



Here we are mainly interested in the continuum limit, and thus the coefficients for v = 0. 
One expects to be equal to the universal one-loop anomalous dimension, = 7q, and 
is the finite part of the one-loop renormalization constant which determines the one-loop 
matching between the SF and the lat schemes, 

*SF& = 4- (5-6) 

We have analysed the series using standard extrapolation techniques |48I49| . In all cases 
we first checked that the coefficients s$ of the logarithm are indeed given by the universal 
one-loop anomalous dimensions. Within the expected numerical precision this is indeed the 
case: for instance in scheme s = 5 and for the O(a) improved action we get 

4/ 7o + = 1.00(1), so ho = 0.997(6). (5.7) 
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S\n QPrlPTTlP 


r + (r — 1 
'0 V c sw — J-j 


r +( r — o) 


r~(r — 1 ) 




1 


-0.2444(2) 


-0.3346(6) 


-0.0635(2) 


-0.0683(5) 


2 


-0.2917(2) 


-0.3818(6) 


-0.1663(1) 


-0.1711(5) 


3 


-0.2360(2) 


-0.3262(5) 


-0.0551(1) 


-0.0599(5) 


4 


-0.3053(2) 


-0.3954(6) 


-0.1935(2) 


-0.1983(5) 


5 


-0.3004(1) 


-0.3905(6) 


-0.1754(1) 


-0.1801(4) 


6 


-0.3043(2) 


-0.3944(6) 


-0.1790(1) 


-0.1837(5) 


7 


-0.2444(2) 


-0.3346(6) 


-0.0635(2) 


-0.0683(6) 


8 


-0.3137(1) 


-0.4038(6) 


-0.2020(1) 


-0.2067(6) 


9 


-0.3088(2) 


-0.3989(6) 


-0.1838(2) 


-0.1885(5) 



Table 1: The one-loop coefficients (|5.5(l of the renormalization constants for the 9 SF schemes, 
for both improved and unimproved Wilson quarks. 



Having passed this check we subtracted the logarithmic divergence using the expected uni- 
versal values Sq = 7 . The finite parts could then be obtained with an accuracy of 4 
significant digits and are given in tabled The precision is generally better for the O(a) im- 
proved data, owing to the fact that O(a) improvement of the action together with tree-level 
improvement of the operators and boundary fields implies the vanishing of the subleading 
coefficients sf. 

As a further check we computed the difference between the results for improved and 
unimproved Wilson quarks. These values must coincide with the ones obtained in pertur- 
bation theory on the infinite lattice. More precisely, the relation between the renormalized 
operators in the lattice minimal subtraction schemes 

(^VA+Av)lat(sw) = (OvA+Av)lat( wilson)'^, w il S on' ( 5 - 8 ) 

can be inferred in two ways, leading, at one-loop order, to the equations 

■y±(l) __ -y±(l) _ -,±(1) -y±(l) _ -,±(1) Q X 

^sw,wilson ^DRED,lat(wilson) DRED,lat(sw) ^SF,lat(wilson) ^SF,lat(sw) ' 

Numerically, we set N = 3 and obtain from section 4 

Wilson = -0.089935(4), X^% lsoa = -0.004845(3). (5.10) 

Indeed, passing via the SF scheme reproduces these numerical values albeit to a lesser pre- 
cision. 

Finally, using the coefficients in tableland combining the results according to eq. (|4.18|) . 
we obtain the two-loop anomalous dimensions in the SF schemes. They are collected in 
table |2J in units of the universal one-loop anomalous dimensions. We observe a large spread 
of numerical values, which already suggests that not all schemes will be well-suited for 
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Scheme 


itht (N = 3) 


7i"/7o ( N = 3 ) 


SF-1 
SF-2 
SF-3 
SF-4 
SF-5 
SF-6 
SF-7 
SF-8 
SF-9 


0.0207(12) + 0.00800(7)iVf 
-0.2394(12) + 0.02377(7)iVf 

0.0669(12) + 0.00520(7)iV f 
-0.3142(12) + 0.02830(7)iV f 

-0.2873(6) + 0.02667(4)iV f 
-0.3087(12) + 0.02797(7)iV f 

0.0207(12) + 0.00800(7)iV f 

-0.3604(6) + 0.03110(4)JV f 
-0.3335(12) + 0.02947(7)iV f 


-0.4668(6) + 0.03890(4)iV f 
-0.1841(3) + 0.02176(2)iV f 
-0.4899(6) + 0.04030(4)TV f 
-0.1093(6) + 0.01723(4)iV f 
-0.1591(3) + 0.02025(2)iV f 
-0.1492(3) + 0.01965(2)iV f 
-0.4668(6) + 0.03890(4)iV f 
-0.0860(3) + 0.01581 (2)JV f 
-0.1360(6) + 0.01885(4)iV f 


DRED 
NDR 
HVDR 


0.093405 - 0.0056290iVf 
-0.011082 + 0.00070367V f 
0.221112 - 0.0133688iV f 


0.045911 - 0.0014072iV f 
0.011082 + 0.0007036iVf 
-0.035357 + 0.0035181iV f 



Table 2: The two-loop anomalous dimensions in units of the corresponding universal one-loop 
coefficients, in various renormalization schemes. 



practical applications. Concerning the equality of the two-loop anomalous dimensions for 
the SF schemes 1 and 7, there is no obvious explanation. In particular we do not see any 
reason why the two schemes should be identical and therefore believe that the equality of 
the anomalous dimensions to one- loop order is an accident. 



6 The step-scaling functions 

Beyond perturbation theory, the running of parameters and renormalization constants is 
traced by computing the corresponding step-scaling functions. For the multiplicatively renor- 
malizable operators (|2,5|) these are denoted by cr ± , and defined in the continuum limit by 



7± 

af(u) = ]imHf(u,a/L), Ef(u,a/L) = 



■ (6.1) 

9 2 {L)=u 



" ~ Z VA+Av . s (g ,a/L) 

Here again, the renormalized quark mass has been set to zero. Similarly, the step-scaling 
function for the running coupling in the SF scheme is defined by 

a{u) =f(2L), u = g 2 (L). (6.2) 

The connection to the renormalization group functions 0, ^ is then given by the two coupled 
equations: 

A„n \y I !„o_ / dg/ _ L (6 .3) 
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L/a 


cltti^ ( L /a) L — i/Cw 


dTfl^ (L 1 — n / Ct? 
uji i vq V / / 1 C S w — U/ w r 


6 


-0.20321867995 


-0.31794582875 


8 


-0.20265948108 


-0.32108637617 


10 


-0.20257791198 


-0.32267336579 


12 


-0.20256208759 


-0.32357398613 


14 


-0.20255806667 


-0.32412959473 


16 


-0.20255683599 


-0.32449503382 


18 


-0.20255642494 


-0.32474770548 


20 


—0.20255629201 


—0. 32492948187 


22 


—0.20255626132 


—0.32506452914 


24 


-0.20255626893 


-0.32516755709 


26 


-0.20255628989 


-0.32524792419 


28 


-0.20255631414 


-0.32531180974 


30 


-0.20255633768 


-0.32536342485 


32 


-0.20255635903 


-0.32540571847 


oo 


-0.20255651209 


-0.32571411742 



Table 3: The one-loop coefficients of the critical mass as obtained from the PCAC condition 
at finite lattice size. For the parameter choices made here, the convergence to the values at 
infinite lattice size is quadratic/cubic in (a/L), for standard/O (a) improved Wilson quarks. 



This implies that the first two coefficients in the perturbative expansion, 

of(u) = 1 + af (1) u + af (2) u 2 + 0(u 3 ), 

are given in terms of bo and 7^,7^ as, 

«f W = 7o ± ln2, 
afW = 7l ± ln2 + 



(In 2)" 



(6.4) 

(6.5) 
(6.6) 



On the lattice the computation of the step scaling functions cr^(n) requires a careful 
extrapolation of lattice approximants T,f(u,a/L), obtained for different lattice sizes L/a 
at fixed values u. The limit is expected to be reached at a rate proportional to a/L, but 
in practice higher order effects may still be important for the accessible lattice sizes. In 
perturbation theory we can now address this question by studying the continuum approach 
of the perturbative coefficients, 

Hf(u, a/L) = 1 + £±(a/L)« u + £±(a/L) (2) u 2 + 0(u 3 ). (6.7) 
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L = 6a 


L = 8a 


L = 10a 


L = 12a 


L = 14a 


L = 16a 




-0.36527 


-0.23780 


-0.15261 


-0.10151 


-0.07018 


-0.05008 


st 


0.02422 


-0.03333 


-0.02885 


-0.01913 


-0.01167 


-0.00653 


st 


-0.48030 


-0.32171 


-0.21396 


-0.14771 


-0.10604 


-0.07863 


si 


0.21731 


0.09522 


0.06117 


0.04699 


0.03880 


0.03320 


si 


0.14610 


0.05147 


0.03183 


0.02601 


0.02308 


0.02099 


si 


0.19678 


0.09254 


0.06316 


0.05017 


0.04211 


0.03631 


si 


-0.36527 


-0.23780 


-0.15261 


-0.10151 


-0.07018 


-0.05008 


si 


0.33235 


0.17913 


0.12252 


0.09319 


0.07466 


0.06176 


si 


0.26114 


0.13538 


0.09317 


0.07221 


0.05894 


0.04955 


si 


0.76256 


0.45542 


0.29780 


0.20823 


0.15308 


0.11690 


S2 


0.31555 


0.20900 


0.14337 


0.10275 


0.07664 


0.05907 


S3 


0.82008 


0.49738 


0.32847 


0.23133 


0.17100 


0.13118 


si 


0.12247 


0.08045 


0.05335 


0.03663 


0.02617 


0.01934 


si 


0.25119 


0.16615 


0.11336 


0.08071 


0.05981 


0.04583 


s e 


0.22927 


0.14607 


0.09736 


0.06810 


0.04975 


0.03765 


si 


0.76256 


0.45542 


0.29780 


0.20823 


0.15308 


0.11690 


s 8 


0.06495 


0.03850 


0.02267 


0.01353 


0.00824 


0.00506 


S9 


0.19367 


0.12419 


0.08269 


0.05761 


0.04189 


0.03155 



Table 4: 0(a) improved Wilson quarks: for finite lattice sizes and SF scheme s = 1, . . . , 9, 5 S 
is defined as the relative deviation of the one-loop step scaling functions from its continuum 
value, cf. eq. (|6,13|) . 



Our one-loop calculation allows the study of 

£±(a/L)« = Zf(2L/a)V - Zt(L/a)M. (6.8) 

This is in principle straightforward, but there is a subtlety related to the determination of the 
zero mass point. With Wilson-type quarks, the chiral limit can only be defined up to cutoff 
effects. Perturbation theory is special in this respect, as go and a can be varied independently 
order by order in the expansion so that a critical mass parameter can be unambiguously 
defined at each order of perturbation theory (cf. (|5.4|) for the one-loop results). However, for 
the perturbative evaluation of cutoff effects to be useful in the analysis of the corresponding 
non-perturbative simulation data, we would like to mimick exactly the procedure used there. 
In particular we take over the definition of the zero mass point: on a lattice of size L/a, with 
T = L and 8 = 0, one determines the bare mass parameter uiq for which the PCAC mass 

, ) = fA(x + a)-f A (x -a) 
m[X0) 4a/ P (x ) (b ' Jj 

vanishes at the midpoint xq = L/2 In perturbation theory, the condition m(L/2) = 
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L = 6a 


L = 8a 


L = 10a 


L = 12a 


L = 14a 


L = 16a 




-1.84058 


-1.59016 


-1.36031 


-1.17808 


-1.03557 


-0.92259 


st 


-1.29860 


-1.15706 


-1.00236 


-0.87333 


-0.77011 


-0.68725 




-1.77678 


-1.55304 


-1.33618 


-1.16118 


-1.02310 


-0.91301 


si 


-1.38914 


-1.20924 


-1.03612 


-0.89692 


-0.78750 


-0.70059 


st 


-1.36068 


-1.19302 


-1.02568 


-0.88964 


-0.78215 


-0.69649 


st 


-1.39431 


-1.21275 


-1.03855 


-0.89868 


-0.78883 


-0.70162 


si 


-1.84058 


-1.59016 


-1.36031 


-1.17808 


-1.03557 


-0.92259 


si 


-1.45294 


-1.24637 


-1.06026 


-0.91382 


-0.79998 


-0.71017 


si 


-1.42448 


-1.23015 


-1.04981 


-0.90654 


-0.79462 


-0.70607 


si 


1.65366 


1.37471 


1.16129 


1.00085 


0.87797 


0.78150 


S2 


1.14358 


0.96017 


0.81540 


0.70456 


0.61875 


0.55094 


S3 


1.62176 


1.35614 


1.14922 


0.99240 


0.87173 


0.77670 


si 


1.23412 


1.01235 


0.84917 


0.72814 


0.63614 


0.56429 


si 


1.17376 


0.97756 


0.82666 


0.71242 


0.62454 


0.55539 


si 


1.19143 


0.98801 


0.83350 


0.71723 


0.62811 


0.55813 


si 


1.65366 


1.37471 


1.16129 


1.00085 


0.87797 


0.78150 


si 


1.26602 


1.03091 


0.86124 


0.73659 


0.64238 


0.56908 


si 


1.20566 


0.99612 


0.83873 


0.72087 


0.63078 


0.56018 



Table 5: The same as table but for unimproved Wilson quarks. The cutoff effects are very 
large, mainly due to the lattice artefacts in the definition of m c on finite lattices (cf. fig. 4). 



then leads to a perturbative series for the critical mass including cutoff effects, 

00 

m c {L/a) = Y J ^ { c ) {L/a)gl\ 

with the low order results 



n=0 



mf\L/a) 
m«(L/o) 



0, 



^(f£\x + a)-f£\x -a)) 



mo=0 



(6.10) 

(6.11) 
(6.12) 



We remark that, in the O(a) improved framework (c sw = 1), the SF correlation function 
/a(xo) in eq. (|6.9|) is supposed to include also the O(a) counterterms proportional to ca and 
Ct- However, at one-loop order and with the chosen parameters, these O(a) counterterms 
vanish identically. The limiting values of m^\L/a) for infinite lattice size L/a are the usual 
one-loop coefficients (|5.4[) . which are reached with a rate proportional to (a/L) 2 and (a/L) 3 
for Wilson and O(a) improved Wilson quarks, respectively. For future reference we have 
collected the values of arni\L/a) for lattice sizes up to L/a = 32 in table 01 
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Having determined the critical quark mass, the cutoff effects in the step-scaling functions 
can be evaluated in a straightforward way. Defining the relative deviation of the one-loop 
coefficients 

5f(a/L) = Et(a/L)V/ofW - 1, (6.13) 

the results for the various SF schemes and both improved and unimproved Wilson quarks are 
given in tables 4 and 5. As can be seen there, cutoff effects in the one-loop coefficient with 
0(a) improved Wilson quarks are typically around the 30-50 percent level at L/a = 6, and 
decrease to a few percent level at L/a = 16. With unimproved Wilson quarks, however, the 
cutoff effects are found to be much larger, typically going from 150 down to 60-80 percent 
at the largest lattice size. We found that most of this dramatic effect is indeed due to the 
usage of m c (L/a) rather than m c (oo). This is illustrated in fig. 4, where the corresponding 
values for 5f and 5% are plotted both for improved and unimproved Wilson quarks. 
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Figure 4: One-loop lattice artefacts 5f and 5^ in the step-scaling function. The data points 
are obtained with the Wilson and SW actions, each for two definitions of the critical mass. 
The dashed lines are obtained as 4-parameter fits to the expected asymptotic behaviour and 
are displayed to guide the eye. 

We conclude that cutoff effects in the step-scaling functions can be quite large, and the 
expected asymptotic dominance of linear lattice artefacts is not yet observed for our data. 
However, we also note that cutoff effects with the O(a) improved action are significantly 
smaller, a fact that is also reflected in the non-perturbative data jSUJ. Moreover, for the 
available lattice sizes 0(a 2 ) effects seem to dominate. We take this as an indication that 
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0(a) improvement of the four-quark operators may be numerically unimportant, at least for 
the step-scaling functions considered here. 

7 Conclusions 

We have introduced a family of finite volume renormalization schemes for the two multi- 
plicatively renormalizable four-quark operators in eq. Q2.5|) . The schemes are based on the 
Schrodinger functional and are defined independently of a particular regularization. By 
matching, at one-loop order of perturbation theory, to commonly used continuum schemes 
(NDR, HVDR, DRED), we could infer the 2- loop anomalous operator dimensions in these 
SF schemes. These results are being used in a corresponding non-perturbative study [SU] 
which completely solves the non-perturbative renormalization problem for these operators 
in quenched QCD. Based on this work, preliminary results for the kaon bag parameter 
Bk have been presented in |12j . Besides the two multiplicatively renormalizable operators 
studied here, the complete basis of parity-odd four-quark operators contains eight further 
operators which form four pairs that mix under renormalization |17j . The study of these 
mixing problems both in perturbation theory and beyond is left for future work. 
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